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ABSTRACT 

We use the photometric redshift method of Chakrabarti & McKee (2008) to infer photometric 
redshifts of submillimeter galaxies with far-IR (FIR) Herschel data ^ obtained as part of the PACS 
Evolutionary Probe (PEP) program. For the sample with spectroscopic redshifts, we demonstrate 
the validity of this method over a large range of redshifts (4 > z > 0.3) and luminosities, finding an 
average accuracy in (1 + z p hot)/(l + z spcc ) of 10 %. Thus, this method is more accurate than other 
FIR photometric redshift methods. This method is different from typical FIR photometric methods in 
deriving redshifts from the light-to-gas mass (L/M) ratio of infrared-bright galaxies inferred from the 
FIR spectral energy distribution (SED), rather than dust temperatures. Once the redshift is derived, 
we can determine physical properties of infrared bright galaxies, including the temperature variation 
within the dust envelope, luminosity, mass, and surface density. We use data from the GOODS-S field 
to calculate the star formation rate density (SFRD) of sub-mm bright sources detected by AzTEC 
and PACS. The AzTEC-PACS sources, which have a threshold 850 /zm flux > 5 mJy, contribute 15% 
of the SFRD from all ULIRGs (L m > 10 12 L Q ), and 3 % of the total SFRD at z ~ 2. 

Subject headings: Galaxies: infrared - Galaxies: submillimeter, photometric redshifts - radiative 
transfer 



1. INTRODUCTION 

Much of the energy emitted by the universe in its in- 
fancy was from dust-enshrouded luminous galaxies. The 
Herschel Space Telescope has opened a new window 
into this epoch, yielding unprecedented sensitivity in the 
far-infrared (Pilbratt et al. 2010), where dusty galax- 
ies emit most of their radiation. The FIR SED allows 
us to probe the physical conditions of dusty sources. 
Chakrabarti & McKee (2005; henceforth CM05) devel- 
oped an analytic means of self-consistently solving the ra- 
diative transfer equation for spherically symmetric, cen- 
trally heated dusty sources. They derived a simple and 
intuitive form for the emergent SED that can be applied 
to infer the physical parameters of dust envelopes from 
the observed FIR SED. Chakrabarti & McKee (2008; 
henceforth by CM08 we refer to the paper and by " CM" 
we refer to CM08's photometric redshift method) ap- 
plied the method of CM05 to fit observed FIR SEDs 
of ULIRGs and SMGs and showed that accurate pho- 
tometric redshifts could be inferred from derivation of 
the light-to-gas mass (L/M) ratio. 

CM08 demonstrated the accuracy of their method with 
a sample of SMGs with FIR data from Kovacs et al. 
(2006). This was the only FIR sample of SMGs available 
at the time, which constituted a sub-sample of the bright 
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SMGs studied by Chapman et al. (2005). The new PACS 
(Poglitsch et al. 2010) observations of SMGs (Magnelli 
et al. 2012) provide an ideal sample for testing the CM 
method on a larger and more diverse sample. SMGs 
are high-redshift galaxies (z > 0.3) classified on the ba- 
sis of their submillimeter flux. "Classical" SMGs have 
-F850 /zm ^ 5 mJy, which corresponds to a luminosity of 
~ 3 x 10 12 L Q at z ~ 2. ULIRGs are generally taken 
to refer to galaxies emitting > 10 12 L© in the infrared 
(Soifer et al. 1984). CM08 had noted that local ULIRGs 
have higher L/M values than their high-redshift cousins, 
a point that we discuss further in §3. The galaxies in 
the Magnelli sample cover four blank fields (GOODS-N, 
GOODS-S, COSMOS, and Lockman Hole) as well as a 
number of lensing clusters. They span a range in lumi- 
nosities from 7 x 1O 13 L > L > 5 x lO n i , and redshifts 
between ~ 0.3 — 5. Many of these sources have spec- 
troscopic redshifts. Thus, this sample is diverse enough 
(Magnelli et al. 2011b) to yield a very robust measure 
of the accuracy of the CM photometric redshift method. 
This is an important test, as photometric redshifts will be 
crucial for analyzing the vast bulk of the observations ex- 
pected from Herschel. Having derived photometric red- 
shifts, we then calculate the star formation rate density 
of sources detected by AzTEC and PACS as a function of 
redshift in the GOODS-S field, which is homogeneously 
covered in the sub-mm down to the classical SMG thresh- 
old (F a5Q /jm > 5 mJy). 

On the basis of IRAS observations, Sanders et al. 
(1988) suggested that local ULIRGs may be produced by 
the merger of gas-rich spirals. The development of hy- 
drodynamical codes that adequately model the collision 
of gas-rich spirals and the subsequent starburst that re- 
sults from a violent merger (Springel et al. 2005) made 
it possible to test this suggestion, which had been an- 
ticipated early on by Toomre & Toomre (1972). The in- 
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frarcd emission of ULIRGs and SMGs during their life cy- 
cles was subsequently calculated using three-dimensional 
self-consistent radiative transfer calculations through the 
time outputs of SPH simulations of merging spirals with 
central AGN (Chakrabarti et al. 2007; Chakrabarti et 
al. 2008; Chakrabarti & Whitney 2009) . These calcula- 
tions reproduced empirically derived correlations, such 
as the correlation between the ratio of the 25 /im to 
60 fim fluxes and energetically active AGN (de Grijp et 
al. 1984; Chakrabarti et al. 2007). Chakrabarti et al. 
(2008) reproduced the clustering of sources in Spitzer 
IRAC color-color plots (Lacy et al. 2004), and explained 
that it was due to the prevalence of the starburst phase 
in the time evolution of these sources. SMGs formed 
during these simulations have diverse properties (Hay- 
ward et al. 2012) and constitute a heterogeneous group. 
However, the simulations analyzed in these papers were 
not cosmological simulations, but rather simulated bi- 
nary mergers of gas-rich systems with central black holes 
that yield high star formation and accretion rates onto 
the central AGN (Springel et al. 2005). Thus, the num- 
ber density of the SMG population as a function of red- 
shift could not be derived. (See however Hayward et 
al. 2011 for a phenomcnological derivation of the num- 
ber density of SMGs in this context) Recently, Dave et 
al. (2010) performed hydrodynamical cosmological sim- 
ulations and identified SMGs as the most rapidly star 
forming systems that match the observationally deter- 
mined stellar mass function of SMGs. In these simula- 
tions, SMGs sit at the centers of large potential wells, 
and accrete gas-rich satellites, but are not typically un- 
dergoing major mergers. However, Magnclli et al. (2010; 
2012) noted that their high star formation rates for the 
bright SMGs (SFR ~ 1000M© yr" 1 ) arc difficult to rec- 
oncile with Dave et al.'s (2010) simulations that produce 
SMGs with star formation rates lower by ~ 3 than what 
is inferred observationally. Thus, models of SMGs can- 
not yet fully account for their contribution to the cosmic 
stellar mass assembly. 

Using our method to derive photometric redshifts from 
the FIR SED can potentially yield a robust measure 
of the number density of SMGs and their contribution 
to the star formation rate density of the universe (and 
thereby their contribution to the cosmic stellar mass 
assembly), without requiring model dependent assump- 
tions. Deriving fundamental parameters of the brightest 
infrared galaxies, such as their luminosity, star forma- 
tion rate, and dust temperature, cannot be accomplished 
without first deriving their redshifts. Unlike optical stud- 
ies (Madau et al. 1998; Hopkins 2004), the derivation of 
star formation rates from the infrared luminosity is not 
affected by extinction, thereby yielding the most robust 
measure of the contribution of dusty galaxies to the star 
formation rate density of the early universe. 

The paper is organized as follows. In §2.1 and §2.2, 
we review very briefly the formalism of CM05 with at- 
tention to the intuitive form of the emergent SED, and 
the redshift inference method presented in CM08. In §3, 
we present the observational sample. In §4, we demon- 
strate the accuracy of the method with respect to the 
Herschel sample with spectroscopic redshifts. In §4.1, 
we apply our method to sources detected by AzTEC and 
PACS in the GOODS-S field and present their star for- 



mation rate density as a function of redshift. We discuss 
future work and present caveats in §5, and conclude in 
§6. 

2. PHOTOMETRIC REDSHIFTS: THEORY 

2.1. SEDs of Dusty Sources 

In CM05, we formulated an analytic solution for the 
FIR SEDs of spherically symmetric, homogeneous dusty 
sources with a central source of radiation. We consid- 
ered envelopes that emit most of their radiation at wave- 
lengths longer than 30 pm, and are sufficiently opaque 
that emission from the dust destruction front does not 
significantly influence the FIR SED. We defined charac- 
teristic parameters, R c h and T c h that are analogous to 
the Rosseland photospheric radius and temperature re- 
spectively, such that 

L = 4nLRl h aTX , (1) 

where a = 5.67 x 10~ 5 erg cm~ 2 s _1 K~ 4 is the Stefan- 
Boltzmann constant. L is a number of order unity that 
allows for better agreement with the numerical solutions, 
and is relevant for very extended atmospheres, which 
have R c ^> I, reflecting the effective increase in emit- 
ting area. We find that 

L = 1.6R" 1 (2) 

is accurate to within ~ 10 % for dust envelopes with 
density scaling as p oc r~ kp with 1 < k p < 2, and 
R c = R c /R ch . Here, R c is the outer radius of the dust 
envelope; our relations apply for R c > 1. 

We assume that the characteristic frequency close to 
which the dust envelope emits most of its radiation, f c h, 
is within the frequency range where the opacity is ap- 
proximately a power law: 

K u = K u a [v/vo) 13 (30 /im < A < 1 mm), (3) 

with j3 = 2 (Wcingartner & Draine 2001). For vq — 
3 THz, corresponding to Ao = 100 /im (fiducial values 
chosen for convenience), we adopt an opacity per unit 
mass of gas of 

>«ioo/mi = 0.54(5 cm 2 g" 1 . (4) 

For 8 = 1, this is twice the value given by the Wein- 
gartner & Draine (2001) (henceforth WD01) model for 
dust in the diffuse interstellar medium since we assume 
that grains in star-forming regions have ice mantles that 
double the FIR opacity. The WD01 opacity is based on 
a gas-to-dust mass ratio of M/Md = 105.1; since the ice 
mantles most likely have a different opacity per unit mass 
than the WD01 grains, we do not attempt to infer the 
dust mass in the sources. Deviations from solar metallic- 
ity, or from the assumed dust model, can be taken into 
account by choosing a different value for S. 

The SED variable R c measures the extent of the dust 
envelope probed by the far-IR SED, which is reflected 
in the shape of the SED: very opaque envelopes have 
R c <~ 1 (close to a blackbody), whereas less dense en- 
velopes can have R c » 1 (very broad SED). We model 
the emergent spectrum by assuming that the emission in 
each frequency channel comes from a shell of thickness 
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Ar m (v) centered at a radius r m (y), with a source func- 
tion (2hv^ h /c 2 ) exp [— hv/kT(f m )] located at an optical 
depth T v (f m ): 



F,. x ( 2 ^) ~ Ku u 3 (k p -l)r 2 rn - h A, „ 



V c 2 



x exp 



hv 



kT(f m ) 



where the optical depth t„ from r to the surface of the 
cloud is 



T v = K v \f kc ' +l 



R„ 



(6) 



Equation [5] (CM05) is the analytic expression that we fit 
to the observed SEDs, where the quantities v = v/v c h, 
r m = r m (v)/R ch = f(R c ,T ch ), and ch = R ch /D A , where 
Da is the angular diameter distance of the source. The 
characteristic frequency v^ = kT^/h, and the quan- 
tity Af m is a function of the SED variables T c h and R c . 
Thus, given three photometric data points we can fit for 
the unknowns T c h, R c and 6 c h in Equation tal This an- 
alytic form for the SED is easily understood as a dust 
envelope with photospheric radius i? c h emitting most 
of its radiation at v^ (close to the peak frequency of 
the SED). Each frequency channel emanates from radius 
T m (y) with thickness Af m ; the high frequency photons 
must also contend with the opacity of the dust enve- 
lope relative to their point of origin and the Wien cutoff. 
CM05 and CM08 have given expressions for the charac- 
teristic emission radii of different frequency channels - 
Rayleigh- Jeans: hv/kT « 1, intermediate: hv ~ kT 
and high-frequency: hv/kT 3> l,r > 1, which comes 
from the outer parts of the envelope, close to the photo- 
sphere, and regions that mediate an effective tug of war 
between the opacity and temperature gradient of the en- 
velope, respectively. The characteristic emission radius 
r m (v) was self consistently determined in CM05 by the 
requirement that the integrated luminosity equal the to- 
tal luminosity. 

Once the SED variables T c h and R c are determined, 
we can solve for the source parameters of the dust enve- 
lope, i..e, the light-to-mass ratio (L/M) and the surface 
density (E) using the following equations, which may be 
inferred from equations 51 and 52 in CM05: 
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M 




(7) 



(8) 



2.2. Inferring the Redshift 

The parameters we infer, L/MS and £<5, depend on 
redshift through the dependence of frequency on redshift. 
We can express the redshift dependence of these param- 
eters in a very simple manner. Since the luminosity of a 
dust envelope satisfies L cx i? 2 T 4 and the inferred mass 
is M oc Ei? 2 , we have: 



L T 4 
MS a E£ 



(9) 



(CM08). The redshift dependence of the surface den- 
sity, £(5, follows from noting that the optical depth at 
the observed frequency must match that at the emitted 
frequency, since r v determines the shape of the SED, 
which is invariant: K(^ obs )(E<J) obs = K(v xest )(ES) IBat , 
so that E obs /E rcst = (^rost/^obs)' 9 = (1 + z) p . Since 



(5) T obs = T rcst /(1 + z), it follows that (CM08) 
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We can now use equation (10), setting (3 = 2 (WD01) 



to infer the redshift of a dusty galaxy from its observed 
value of L/MS, which we have shown can be derived from 
the FIR SED analytically (CM08): 



1 + Zinf = 



(L/MS) 



"I 1/6 



(V^)obs 



(11) 



where (L/MS) is the typical value in the 
Herschel sample; since L/MS can range over an 
order of magnitude, we use the geometric mean. Note 
that L is the luminosity derived from fitting to the 
mm-FIR SED. A simple way to understand why our 
method yields a more accurate redshift than estimating 
it from the temperature (or peak of the SED) is as 
follows. Our inference of the redshift is proportional 
to L/M 1 /®. Alternately, inferring the redshift from 
fitting a blackbody function to the SED and estimating 
a temperature from the shift of the peak of the SED 
is proportional to the temperature itself (i.e., to the 
template rest wavelength temperature assumed). It 
follows from equation that L/M oc T^R- 1 (for 

k p = 2 and (5 = 2). The parameters (T c h and R c ) 
are correlated inversely for a given luminosity (from 
equation 1), leading to a smaller range in L/M 1 / 6 than 
in the temperature itself, and therefore a more accurate 
redshift. This is the case when R c can be independently 
determined, which can be done when there is sufficient 
information on the shape of the SED, as we discuss in 
§4. 

3. OBSERVATIONS 

In this study, we used deep PACS 70, 100 and 160 jum 
observations provided by the Herschel Space Observa- 
tory as part of the PACS Evolutionary Probe (PE#]) 
guaranteed time key program (Lutz et al. 2011). Specif- 
ically, we used here the observations of the Great Ob- 
servatories Origins Deep Survey-North (GOODS-N) and 
South (GOODS-S) fields, the Lockman Hole (LH) field, 
the Cosmological evolution survey (COSMOS) field and 
the lensed fields Abell 2218, Abell 1835, Abell 2219, Abell 
2390, MS1054, CL0024 and MS045. The sample with 
spectroscopic redshifts is taken from and fully described 
in Magnelli et al. (2012). 

The first aim of this paper is to test the CM FIR pho- 
tometric redshift method on a large and more diverse 
sample than studied by CM08. To calibrate and test this 
photometric redshift method we use a sample of SMGs 
detected by PACS and with robust spectroscopic esti- 
mates from Magnelli et al. (2012). Here, we restrict our 

6 http: / /www.mpe.mpg.de/ir/Research/PEP 
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study to the sub-sample of PACS detected SMGs and 
we do not use the SPIRE flux densities in the compu- 
tation of our photometric redshifts. Our final sample 
contains 80 PACS detected SMGs, which have sub-mm 
(850 /im, 870 /im) or (and) mm (1.1 mm) data (Magnelli 
et al. 2012). Of these sources, about one-half have ro- 
bust spectroscopic redshift estimates. As in CM08, we 
restrict our SED fitting to sources that have at least three 
photometric data points with S/N > 3. 

As a second step we apply the CM photometric method 
to a PACS-detected SMG sample of the GOODS-S field 
(i.e. without requiring any spectroscopic redshift esti- 
mate) and derive their SFRD. This SMG sample corre- 
sponds to the AzTEC (1.1 mm) sample of the GOODS-S 
field presented in Scott et al. (2010) and which con- 
tains 41 sources detected with S/N> 3.5. This AzTEC 
sample is cross-matched with our PACS catalogue us- 
ing a matching radius of 9". We find that 25/41 AzTEC 
sources are detected in our PACS catalogue. We note as- 
suming a 5g5o / jm/<S'iioo Mm ratio of 2.1 (Scott et al. 2008; 
Austermann et al. 2010), the depth of the AzTEC sur- 
vey converts into a 850 /im flux of 4.2 mJy Thus, us- 
ing our final sample of 25 galaxies, we derive the SFRD 
of AzTEC and PACS detected sources (i.e., down to 
~ 1.2 mJy and ~ 2.4 mJy at 100 and 160 /im respec- 
tively) sources that have 850 /tm fluxes typical of " clas- 
sical" SMGs (i*85u urn > 5 mJy). Our results on the 
SFRD may be uncertain in part due to the PACS non- 
detections of some of the SMGs in the GOODS-S field, as 
well as the relatively small number of sources. Therefore, 
we emphasize that we derive here the SFRD of a sample 
of infrared bright galaxies that have extrapolated 850 /im 
fluxes typical of classical SMGs, specifically the AzTEC- 
PACS sources in the GOODS-S field, rather than the 
SFRD of the SMG population. 

4. PHOTOMETRIC REDSHIFTS FROM A SAMPLE OF 
HERSCHEL SOURCES 
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Fig. 1. — Accuracy of redshift determination using the CM pho- 
tometric redshift method; 2; n f is the photometric redshift we de- 
rive and z is the spectroscopic redshift. Yellow filled circles mark 
sources with long wavelength data that sample close to the break 
frequency. 

We first examine the sub-set with spectroscopic red- 
shifts to verify the accuracy of the CM method. Some 



of the sources in this sample are lacking millimeter data, 
i.e., data points that are close to the break frequency of 
the SED, where the character of the emission switches 
from being dominated by hotter photons coming from 
close to the photosphere to colder photons that come 
from the outer regions of the envelope. The break fre- 
quency of the SED is given by (CM08): 

iW = 6.35 x 10 9 (^) ' -^Hz . (12) 
V R c J t + z 

For average values of T c h ~ 100 K and R c ~ 300, galaxies 
at z ~ 2 would have a break wavelength of ~ 2 mm. We 
have determined (by successively removing long wave- 
length data points of local ULIRGs and fitting the SEDs) 
that if the lowest frequency data point is more than a 
factor of two higher than the break frequency, then the 
shape of the entire SED (including the transition to in- 
termediate frequencies) can not be clearly determined. 
In such cases, we adopt a median value of R c (R c = 300) 
that is determined by fitting to the sample that has both 
spectroscopic redshifts and long wavelength millimeter 
data, and at least three data points with S/N > 3 (with 
the latter being a requirement for all the sources we fit). 
We note that the variation in the dust temperature when 
the SED is fit with a blackbody function is slightly larger 
than the variation of the characteristic temperature when 
the SED is fit with our model. In some cases, as shown 
in Figure [TJ we do have long wavelength data points that 
sample the break frequency of the SED. This allows us to 
improve the accuracy of our photometric redshift method 
(as it would allow us to independently fit for not only T c h 
but also R c ), and also provides a better determination of 
the mass of the envelope. 

The average L/M of the sample with spectroscopic red- 
shifts is 30 Lq /Mq , which is comparable to the average 
L/M determined for SMGs studied in CM08. This is the 
template L/M value that we employ in deriving photo- 
metric redshifts here. CM08 had noted that the L/M 
they derived for SMGs was a factor of two lower than 
local ULIRGs. The somewhat lower values of L/M for 
SMGs relative to local ULIRGs found by CM08 is re- 
flected in the lower dust temperatures found in high red- 
shift infrared bright galaxies (Symeonidis et al. 2011). 
Figure [l] depicts the accuracy of our method for sources 
in this sample with spectroscopic redshifts. The yellow 
filled circles mark the sources that have sufficient long- 
wavelength data, i.e., v m i n < 2 x J-Wak, to determine 
R c independently. The average accuracy for the entire 
sample is 1.11, and the standard deviation is 0.13. The 
corresponding numbers when the redshift is inferred by 
fitting a single temperature blackbody function to the 
SED are 1.2 and 0.2 respectively. Here we have used a 
temperature of 60 K, finding that it yielded the maximal 
redshift accuracy when fitting with a blackbody func- 
tion. Using a modified blackbody fit (with (3 = 1.5 as 
recommended by Gordon et al. 2010 and T=50K) yields 
similar accuracies as for a single temperature blackbody. 
The sample with sufficient long-wave data has an aver- 
age accuracy of 1.07 and standard deviation of 0.075. 
The higher accuracy of the CM method is due to fit- 
ting for two parameters (R c and T c h) to account for the 
influence of a range of temperatures on the FIR SED, 



along with the peak frequency of the SED. We can fit 
for both parameters independently when we have data 
for the long wavelength part of the SED, sampling close 
to the break frequency. For the sample that does not 
have sufficient long-wavelength data, we find a smaller 
dispersion of characteristic temperatures when we fit our 
functional form to the data than when we fit a single 
temperature blackbody, as our functional form is a bet- 
ter approximation to the real FIR SED. We attribute 
the slightly higher accuracy of the CM method for the 
sample with restricted data to the smaller dispersion in 
temperatures when we fit the data with our functional 
form. 

It is important to contrast our method with other pho- 
tometric redshift methods that utilize the FIR SED. We 
have already discussed above the accuracy in the infer- 
ence of the redshift by fitting a single temperature black- 
body. Another commonly used FIR redshift method em- 
ploys the radio-sub-mm spectral index (Carilli & Yun 
1999). This technique provides a gross indication of the 
redshift by using the variation of the radio-sub-mm spec- 
tral index as a function of redshift. The uncertainty in 
redshift for this method depends on the scatter of the 
spectral index for a given polynomial fit to the data. 
The uncertainty in redshift also depends on how well a 
given polynomial fits the data. The ±ltr variation in the 
spectral index for the adopted polynomial fit in Carilli & 
Yun (1999) is ~ 1, with 30% of the galaxies in their sam- 
ple falling outside this range. Thus, this method is less 
accurate than ours. Aretxaga et al. (2005) sample mock 
galaxy catalogs with an assumed evolutionary history to 
derive photometric redshifts from radio to FIR colors, 
finding a redshift accuracy of Sz ~ 0.3 (rms dispersion 
about the line z p hot = z spcc ). In terms of our metric 
(1 + Zjnf)/(1 + -z spoc ), (where we take Zj n f to be the most 
probable photometric redshift) the average accuracy of 
the sample considered in Aretxaga et al. (2005) is 1.213 
(the average accuracy of our sample is 1.12). Kovacs et 
al. (2006) suggest the use of the FIR-radio correlation 
as a photometric redshift estimator, in a manner sim- 
ilar to Carilli & Yun (1999). Roseboom et al. (2012) 
use radio, mm and SPIRE data to find comparable av- 
erage accuracy (although the standard deviation in their 
estimates is greater). In summary, these photometric 
redshift methods (whether estimating the redshift from 
Monte Carlo sampling of mock catalogs, the radio-sub- 
mm spectral index, or from the dust temperature) yield 
comparable or generally poorer accuracy relative to our 
method. 

We conclude therefore that our method yields reason- 
ably accurate photometric redshifts for infrared bright 
galaxies out to z ~ 5. Maximal accuracy of our method 
can be realized with the addition of long wavelength mil- 
limeter data sampling the break frequency, so that T c h 
and R c are independently determined. The addition of 
near-IR data (Daddi et al. 2009) does improve the ac- 
curacy of photometric redshifts beyond what we are able 
to do here with FIR data alone. We should note though 
that our method does not require radio data, as is used 
in the estimates by Aretxaga et al. (2005), Carilli & Yun 
(1999), and Roseboom et al. (2012), or an optical identi- 
fication, as is used by Daddi et al. (2009). This may be 
an advantage compared to other methods (those that use 



near-IR data) that achieve similar or better accuracy. 

4.1. Cosmic Star Formation History of AzTEC-PACS 
sources in GOODS-S 
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Fig. 2. — Star formation rate density of AzTEC-PACS sources in 
GOODS-S using Herschel data (black symbols). The yellow and 
green points are the SFRD's derived by Magnelli et al. (2011b) 
using 24 /im and 70 (im data for ULIRGs and all galaxies respec- 
tively in GOODS-S and GOODS-N fields. The shading indicates 
the spread in derived values from Magnelli et al. (20011b) 



Now that we have derived the photometric redshifts, 
we use them to calculate the star formation rate den- 
sity (SFRD) for the PACS selected SMC sample in the 
GOODS-S field. We focus on the GOODS-S field here 
as that field is homogeneously surveyed down to 5 mJy 
at 850 /jm (thereby selecting "classical" SMGs), as dis- 
cussed Jpl The other fields (GOODS-N, Lockman Hole, 
and COSMOS) are surveyed at different depths over 
parts of the field in the mm and sub-mm, and as such 
yield a more heterogeneous sample of galaxies than the 
GOODS-S field observations. 

We convert our infrared luminosities to star formation 
rates using the relation presented in Kennicutt (1998), 
modified (by a factor of ~ 2; Portinari et al. 2004) to 
use the Chabrier IMF (Chabrier 2003), i.e., we take 



SFR = 



1 x 10 10 L Q 



(13) 



The average star formation rate for galaxies with infrared 
luminosity in excess of 3 x 10 12 _Lq is 700 M Q /yr for all 
of the fields we have studied here, including GOODS-S 
(this is consistent with the observed IRLF at z ~ 2, e.g., 
Wardlow et al. 2010). This is a factor of ~ 3 in ex- 
cess of the average values found by Dave et al. (2010), 
and prior to Herschel observations the total infrared lu- 
minosity and hence the star formation rate of SMGs was 
indeed uncertain. However, recent Herschel observations 
which probe the rest-frame FIR peak of SMGs (Magnelli 
et al. 2012) do support the higher star formation rates of 
major merger models (Chakrabarti et al. 2008; Hayward 
et al. 2011). 

Figure 2 shows our derived SFRD (black symbols) 
over-plotted with the determination by Magnelli et al. 
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(2011b) using 24 /xm and 70 /zm data from the GOODS- 
S and GOODS-N fields. The green and yellow symbols 
depict the SFRD for all galaxies at a given epoch and 
all ULIRGs respectively. Note that the sample in Mag- 
nelli et al. (2011b) is a more diverse sample than the one 
we analyze here, including LIRGs and normal galaxies 
as well as ULIRGs. We have derived the SFRD shown 
here (and corresponding infrared luminosity density) us- 
ing the standard 1/Vm ax method (Schmidt 1968). The 
luminosity density is given by: 

4>(z j )Az j = £ (14) 

Here, V c>i = min[U c (z maXii ), V c (Zj)] - V c (Zj-t), with 
V c being the co-moving volume; i labels sources and 
j redshift bins, so that z max ,i is the maximum redshift 
at which source i would have been detected, Azj — 
Zj ~ Zj—i, Li is the luminosity of source i, and Zj = 
0.5(zj_i + Zj). The star formation density in a redshift 
bin is then simply the luminosity density scaled by the 
conversion factor given by Kennicutt (1998), using the 
Chabrier IMF correction. 

We find that the sub-mm bright PACs sources con- 
tribute 15 % to the SFRD at z ~ 2 relative to ULIRGs 
(as determined in Magnelli et al. 2011b), and 3 % rela- 
tive to the total SFRD at that epoch. These percentages 
are relative to the most probable value for the SFRD 
determined by Magnelli ct al. (2011b), with the spread 
in derived values indicated in Figure [2l For comparison, 
Magnelli et al. (2011b) found that the ULIRG popula- 
tion contributes 17 % to the total SFRD at z ~ 2. As dis- 
cussed in £j3j the sources we consider here have a 850 /im 
flux typical of " classical" SMGs (F 850 ^ m > 5 m Jy) . 
However, given the relatively small number of sources 
(~ 30), we interpret our findings as yielding the SFRD 
of AzTEC-PACS sources in GOODS-S, rather than that 
of the SMG population. It is worthwhile noting that 
there are systematic differences in the determination of 
the SFRD of SMGs using different surveys. For instance, 
radio-detected LESS SMGs in the Extended Chandra 
Deep Field South (ECDFS) as studied by Wardlow et 
al. (2010) have a different redshift evolution and over- 
all contribution to the SFRD than the sample studied 
by Chapman et al. (2005). Our results are comparable 
with the SFRD of mm-selected sources in the Lockman 
Hole and GOODS-N fields calculated by Roseboom et al. 
(2012). Our determination of the SFRD of the AzTEC- 
PACS sources is a lower limit to the ULIRG population 
which are a less extreme version of infrared bright galax- 
ies and increasingly common at high redshifts (Wuyts et 
al. 2011). The shape of the SFRD of the AzTEC-PACS 
sources does not show a decline out to z ~ 4. 

5. DISCUSSION 

There are a few caveats worth mentioning, relating 
to the use of our method and astrophysical effects. If 
lensed galaxies are present in our sample, they would ar- 
tificially boost the number counts of SMGs. Negrello et 
al. (2007) estimate the probability of lensing as a func- 
tion of the sub-mm flux of a galaxy, which we can use 
to ascertain the likelihood of galaxies in our sample be- 
ing lensed. They find a significant incidence of lensing 
(> 50%) for galaxies with 850 /im fluxes in excess of 



40 mJy, but all our sources are well below this thresh- 
old. For known lensed galaxies, we have checked that our 
method works equally well on the lower luminosity SMGs 
in the A2218 lensed cluster. We conclude therefore that 
lensing is not a serious source of concern for this sam- 
ple, i.e., the probability of galaxies being lensed in the 
fields we have studied is very low given their 850 fjxa 
flux. Moreover, we have checked that for known lensed 
galaxies in the A2218 lensing cluster, our method yields 
robust photometric redshifts. 

Our simple model is designed to yield through SED- 
fitting the large-scale parameters of infrared bright galax- 
ies, namely the mass, size, and luminosity once the red- 
shift is determined using the CM photometric redshift 
method. We do not treat the effects on the SED of a 
clumpy dust envelope or distributed sources of luminos- 
ity. Although we do not treat the morphology of infrared 
bright galaxies in detail, our derivation of the mass, size, 
luminosity and density profile is in good agreement with 
observational determinations (CM08). The advantage of 
our model is that it incorporates a self-consistent range 
of dust temperatures in an analytic expression for the 
emitted luminosity, which we have shown can be applied 
to infrared bright galaxies to robustly derive redshifts 
and large-scale source parameters out to z ~ 4. 

6. CONCLUSION 

• We have applied the photometric redshift method 
of CM08 to Herschel data to demonstrate the accuracy 
of our method for a large and diverse sample of SMGs. 
We find accuracies of ~ 10% relative to spectroscopic 
redshifts, i.e., in (1 + z in f)/(l + z), out to z ~ 4. 

• We also give the average values for the L/M ratios of 
SMGs in these samples, which is lower than that for local 
ULIRGs by a factor of 2. Our average star formation rate 
is 700 Mq /yr for galaxies with luminosities in excess of 
3xl0 12 L Q . 

• We estimate the star formation rate density of sub- 
mm bright PACS sources in the GOODS-S field. These 
sources have an extrapolated 850 /im flux typical of clas- 
sical SMGs, ^850 ^m ^ 5 mJy (we extrapolated the sub- 
mm flux from AzTEC observations as discussed in q3j). 
Our derivation of the SFRD is a lower limit, particularly 
at low redshifts [z < 1) where normal spirals (and secular 
processes or minor tidal interactions) drive the star for- 
mation history of the universe. These PACS sources con- 
tribute 15 % to the SFRD relative to ULIRGs at z ~ 2, 
and 3 % relative to the total SFRD at that epoch that is 
produced by all galaxies. We find that there is no decline 
in the shape of the SFRD of the sub-mm bright PACS 
sources in the GOODS-S field out to z ~ 4. 
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